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O ; Abstract 

We discuss new possibilities for Off-Diagonal Long Range Order (ODLRO) 
K^ in spin chains involving operators which add or delete sites from the chain. 

j~^ ' For the Heisenberg and Inverse Square Exchange models we give strong nu- 

o' 

merical evidence for the hidden ODLRO conjectured by Anderson [|[|. We 
find a similar ODLRO for the XY model (or equivalently for free fermions 
C^ ' in one spatial dimension) which we can demonstrate rigorously, as well as 

numerically. A connection to the singlet pair correlations in one dimensional 
models of interacting electrons is made and briefly discussed. 
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Typeset using REVT^ 



In 1991 one of us (P.W.A) conjectured, based on RVB ideas for the one-dimensional 
Hubbard model, that there should be a non-zero overlap between the groundstate of the one 
dimensional Heisenberg model on a chain of A^ sites and the state obtained by inserting a 
pair of nearest neighbor spins in a singlet configuration into the ground state of the A^ — 2 site 
Heisenberg chains |]I],0]. Both of these models were therefore expected to have a hidden form 
of ODLRO for an operator which not only involved sampling the state of the system (in this 
case checking that a given pair of spins was in a singlet) but also changing the Hilbert space 
of the system in an essential way (adding two sites), reminiscent of the Girvin-MacDonald- 
Read order parameter of the Fractional Quantum Hall Effect (FQHE) [||, despite the fact 
that that the Heisenberg model has gapless excitations. The conjecture that the overlap 
should be finite has not been previously investigated (but see Q] where evidence for the 
resulting ODLRO was found); therefore we numerically tested the original conjecture that 
the overlap for the ISE and Heisenberg models between the N site groundstate and the 
N — 2 site groundstate with a local singlet inserted should be non-zero. 

These overlaps as a function of A^ are shown in Fig. |I] together with fits to the results 
of the form 0.817 + 0.778A^~2 fQ^ the ISE model and 0.820 + .7A0N~^ for the Heisenberg 
model. The results strongly suggest that both overlaps remains finite in the limit as the 
system size goes to infinity. Note that since the phases of the A^ and N — 2 site wavefunctions 
may be chosen independently the phase of the overlaps is meaningless and further, since the 
groundstate momenta of the A^ and N — 2 site groundstate wavefunctions differ by vr, the 
overlap is multiplied by minus one if the location of the singlet pair is shifted by one site. 
We have chosen the overlap real and positive for convenience. 

As Fig. U shows, the overlaps in the ISE and Heisenberg models are not only finite in the 
N —>■ oo limit, but also surprisingly close to each other, despite the fact that the ranges of the 
interaction in both models are quite different. We now present an analytical calculation in 
the ISE model that gives some understanding of origin of this finite number. The groundstate 
wavefunction of the ISE model in a basis of local spins {|o"i ■ ■ ■ <Jm)}, ctj = ±i is given by 0: 
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For the A^ site ISE groundstate M = N and {zi} = Cn = {e^^j^^i, while for the A^ - 2 
site groundstate M = iV — 2 and {zi} = Cn~2 = {^'^~^}n=i To compute the ISE overlap 
we add a^^i, ctat, sitting in a singlet, to the N — 2 site groundstate. This overlap is hard 
to calculate since the Zj's from both sets are not commensurate with each other. However, 
if we slightly deform the set CAr_2 to be {e^^j^jf and leave (Jn-i, cfn in a singlet then we 
obtain a new state \l/^, that can be recognized as a localized 2-spinon state f^J^l- Spinous 
are the elementary excitations of the ISE model with semionic statistics ||^. Although this 
localized spinon state is not an energy eigenstate it consists of an admixture of eigenstates 
that contain only spinous (the groundstate on A^ sites) or 2 spinous. Here the localized 
spinous sit at sites A^ — 1 and A^ in a singlet. In general we could have put them at sites 
a, P by deforming Cn~2 into {e^^\i = 1, . . . , N}/{e^^ , e^^}. 

In the basis of states with M overturned spins with respect to the ferromagnetic state 
labeled by their positions along the chain: rii, . . . , um, ^at and \l/^ are given by: 
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We now calculate the overlap between \1/^ and \I'^ for arbitrary separation a between 
the spinous (because of translational invariance we can fix f3 to be 0). We have for the 
overlap: 

(^0|^2(a)) _ (^0|^2(a)) f {^^{a)\^^{a))Y'^ 



First we will determine (\I'°|\I'2(a)). Setting M = f - 1: 
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In eq. (|^) we shifted all the rii by a in the second term to bring both terms in the same 
form. 

Since \E'° is a singlet we know that if the two spinons in \E'^ were in the 8^ = triplet 
state the overlap should be zero, i.e. (\E'°|\E'^ (a), triplet) oc a{a) (1 + (— )") = for all a. 
Therefore a (a) vanishes for all even a. At the same time we see from eq. @ that a {a) is 
a polynomial in cos (t^)- This is easily checked by expanding the sin f "•y^ Trj and noting 
that a{a) is even in a. We conclude immediately: 

N_ 

a(a)ocj2(cos(^)-cos(^)) 
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Now we turn our attention to the second piece. Only its numerator is a dependent: 

2(^-1)' E (^'(ni,...,nM|0,a))' 

{ni,...,njvf} 

But this sum can be recognized as the {Sz{a)Sz{0)) static correlation function in the ISE 
model 0, or — after an (exact) conversion of the sums to integrals — as the one-particle 
density matrix in the Calogero-Sutherland model at half filling f^ . The result for large A^ is 



that is proportional to ''■^"^ , where Si(x) is the sine- integral function. The entire expression 
for the overlap, with normalization computed by considering a = 0, then becomes: 

2 sin (^) / vra 
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Thus for a nearest neighbor the overlap is ^\/g^ — 0.82917, which is within 1.5% of the 
Heisenberg and ISE groundstate singlet insertion overlaps. 

Many of the properties of the two spinon matrix element can be understood in terms of 
the Girvin-MacDonald-Read P] order parameter for the bosonic i^ = \ Laughlin [^ state. 
In the spinon language, the insertion of up-spin spinous acts on the down-spin wavefunction 
in that same way as the quasi-hole operator of the i^ = \ state, while the insertion of a 
down-spin spinon acts analogously to the insertion of a quasi-hole and a hard-core boson. 
Consequently, the two spinon insertion for a = clearly takes the bosonic ^ = -^ state 
with Y — 1 particles, periodic on a ring of size N — 2, to the i^ = ■^ state with y particles 
periodic on a ring of size N, and is essentially the Girvin-MacDonald-Read order parameter 
for the i^ = I state and the overlap associated with it must be identically equal to one. 
Further, for a 7^ 0, the two spinon insertion is equivalent to a sum of two operators: a 
quasi-hole being inserted a finite distance from a hard core boson and another quasi-hole. 
Such an object should still have an expectation value given essentially by the quasihole 
propagator. This would decay exponentially in the bulk of a z/ = ^ state, but near the edge 
would decay only as a" 2, exactly as the two spinon matrix element does (see Eq. (^) and 
Fig. 1^). The vanishing of the matrix element for any even a and its alternation for odd a 
simply reflect the lattice structure of our model, the phase of the quasi-hole propagator and 
the fact that the spinon operator is a superposition of two operators so that the result is 
real, rather than complex. The reason why the two spinon overlap agrees so well with the 
actual singlet insertion overlap is unclear, however, our numerical Monte Carlo results for 
the singlet insertion display the same power law decay as a function of the spin separation 
(see figure 0). The insertion overlap also vanishes for separations which are even multiples 
of the lattice spacing and alternation for separations which are odd multiples of the lattice 



spacing, so it would appear that for the Heisenberg and ISE models the spin chain ODLRO 
is essentially that same as the ODLRO of the i^ = ^ Laughlin state. 

Note that, since the two spinon insertion operator generates states with either zero or two 
spinons, the zero spinon state being the ground state, its action is analogous to the action 
of a pair of fermionic creation and annihilation operators on the fermionic groundstate, 
generating states which contain one particle and one hole or else no particles and holes, 
i.e. the groundstate. The matrix element to the ground state as a function of the space- 
time separation of the fermionic creation and annihilation operators is the amplitude for 
the fermion to propagate to the space-time location where it can be annihilated and is by 
definition the fermion propagator. The overlap we calculate is, by analogy, proportional to 
the equal time spinon propagator, i.e. the amplitude for either of the spinons to propagate 
to the location of the other, at which point the two can annihilate each other, provided they 
are in a singlet configuration. 

In an effort to better understand our results for the Heisenberg and ISE models, we 
examined the overlap for the same insertion operator acting on the N — 2 site XY model 
groundstate with the A^ site XY model groundstate. The strikingly similar result, here 
extended to a larger system size than was possible for the Heisenberg and ISE models, is 
shown in Fig. |l]. Since the XY model can be mapped onto spinless fermions, the "order 
parameter" can be recast in that language and a more detailed study made. In that lan- 
guage the non-zero overlap is between the groundstate for N/2 spinless fermions, with either 
periodic or antiperiodic boundary conditions, on A^ sites (with positive hopping integral so 
that their momenta in the groundstate are centered around vr) and the state obtained by 
adding two sites and one fermion (in a superposition of being on the two added sites with 
a relative minus sign for the two different sites) to the groundstate for the N/2 — 1 spinless 
fermions on A^ — 2 sites — with the opposite boundary conditions from the A^ site case. The 
change of boundary conditions is a non-local operation in the fermion language, however in 
the spin language it arises from a local operation, i.e. the insertion of two additional spins 
in a singlet configuration. 
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Since the wavefunctions for the two fermion states to be overlapped are those of free 
particles, considerable progress can be made in its computation. In particular, it can be 
bounded rigorously from below. The calculation is straightforward, using simple proper- 
ties of polynomials which roots lie on the unit circle, and the fact that the corresponding 
momenta for the N — 2 and A^ site models differ only by 0{1/N'^) near the Fermi surface. 
We will present the proof elsewhere. Here we state only the result that the lower bound 
obtained was e~2 + 0{N~^). 

The non-zero overlap for the XY model is particularly striking since we know that the 
overlap between the A^ site, y particle groundstate and the state created by inserting a 
single, localized electron into the A^ site, y — 1 particle groundstate would vanish like iV~2, 
while the overlap between the A^ site, y particle groundstates between models with periodic 
and antiperiodic boundary conditions vanishes like A^~^. 

The vanishing of the overlap for the insertion of a single, localized fermion implies that 
no analog of the two spinon calculation exists for the XY model, since the only way to make 
the momenta similar in the two overlapped states is to make ^ — 1 of them identical is which 
case the overlap will vanish like A^~2. 

It may appear that, since the "order parameter" we have defined involves spin chains with 
different numbers of sites, there will be no physical consequences to this form of ODLRO 
and in fact there is, for example, no additional groundstate degeneracy associated with the 
existence of this kind of order. However, order of the type we find can have important 
consequences for more general models than spin chains. For example, we will now show 
that the ODLRO of the Heisenberg model is responsible for the leading contribution to the 
singlet pair susceptibility of the one dimensional Hubbard model and is thus not without 
potentially important physical consequences. 

First, we make the connection between the ODLRO found for the Heisenberg model and 
the equal time, singlet pairing correlation function of the one dimensional Hubbard model 
in the limit as, U ^ oo. In that limit the groundstate wavefunction of the Hubbard model 
is given by a product of spin and charge wavefunctions, the latter being given by a spinless 
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ferinion determinant and the former by a Bethe Ansatz wavef unction for the "squeezed" 
Heisenberg model, i.e. a Heisenberg model defined only on those sites occupied by the 



spinless fermions [|TT| . For periodic boundary conditions and the number of electrons equal 
to 4N + 2, A^ an integer, the ground state wavefunction for the Heisenberg model should 
be used. The equal time singlet pair correlation function is given by the overlap of the 
groundstate with the state obtained from the groundstate by removing a nearest-neighbor, 
singlet pair of electrons at sites which we can take to be and 1 and then inserting a earest- 
neighbor, singlet pair of electrons at sites j and j + 1. Due to the hidden ODLRO of the 
Heisenberg model, the spin wavefunction overlap for a fixed charge-configuration has a piece 
which for large separation, j, is just given by a constant times (— 1)"^, where Uj is the number 
of spinless fermions found on sites between 1 and j in that charge configuration. This leads 
to a contribution to the correlation function which is given by the expectation value in the 
spinless fermion groundstate of \E'^(j + l)\E'^(j)^(l)^(0)(— l)^j>'>i , similar to the 

correlation function for the alternating, spin-spin correlation function studied by Sorella, 
et al. [|I2|. The leading asymptotic behavior of this expectation value can be computed 



straightforwardly from Abelian bosonization and is given by A co^{kpx)x~2^ where kp is 
the kp of the spinless fermions and A is a cutoff dependent constant. This agrees with the 
predictions of the Luttinger liquid description of the Hubbard model [|l^,0 , provided that 



we take the charge sector rediagonalization parameter, Kp, to be 1/2 and remember that the 
kp of the spinless fermions is twice that of electrons with spin at the same filling fraction. 

In the Hubbard model, this contribution to the singlet pair correlation functions arises be- 
cause the bosonized form of '?/'t(j)'?/'j,(j + 1) contains an operator proportional to exp (i0_R,p), 
involving only charge degrees of freedom. [|13| This operator is present because the operator 



product expansion for exp (|6i?,cr(a;) j exp ( — |0_R,(T(a;') j contains the identity times a coeffi- 
cient asymptotically proportional to (a; — x')~2. The decay of this coefficient with separation 
implies that, if the electrons are inserted m sites apart where m ^ 1 but j ^ m, the singlet 
pair correlations decay with j in the same way, but there is a multiplication of the prefactor, 
A, by a factor of m~ 2 arising from spin degrees of freedom |jT^. This is in in agreement with 
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the decay we find for tlie singlet and two spinon insertions in tlie ISE and Heisenberg models, 
and suggests the identification of the insertion of an up-spin spinon into the ISE model with 
the action of the operator Omsert = V exp i^QR^aiJ o)) + {—iY exp {i2,'^L,a{j ct)), where 
a is the lattice spacing and j is the number of sites to the left of the insertion site in the 
original chain. Since this operator is a semion, it is natural to identify it as the spinon cre- 
ation operator; an identification compatible with the observation of |jl5| that the generalized 
commutation relations of the Fourier modes of this operator provide a natural realization of 
the Yangian. 

The identification is also compatible with the alternation with odd separation and van- 
ishing for even separation that we find for the singlet and two spinon insertions. Both are 
in agreement with the singlet pair correlations of the Hubbard model as computed with 
Abelian bosonization. 

The connection to Abelian bosonization should be extendable to generalizations of the 
Hubbard model which include spin-dependent interactions so that the Luttinger liquid re- 
diagonalization parameter of the spin sector of the Hubbard model, k^, (for a definition of 
Ko- see |]T6[) is renormalized from 1. At half filling the low energy sector of the model then 
becomes a general XX Z model. This would change the exponent for the decay of the singlet 
insertion overlap from | to an exponent, rj = ^{na + i^a^)- For XX Z models this is a simple 
function of the anisotropy of the model 0]. In general: k^ = 1 — ^ cos~^ (j^)- This is 
equal to | at the XY point exactly as we observe. Based on the connection between the two 
spinon overlap and the singlet insertion for the ISE model, where spinous are noninteracting, 
we argued that the decay exponent was just the exponent of the spinon-spinon propagator. 
In this language, the continuously varying exponent of the overlap decay in general XX Z 
chains is just the exponent of the spinon-spinon propagator in models of interacting spinous. 

We note that that the alternating piece in the singlet-pair correlation function is the 
slowest decaying piece of that correlation function for the Hubbard model and may have 
important consequences for superconducting correlations in models based on the one dimen- 
sional Hubbard model. In particular, if one considers an array of Hubbard chains coupled 



with operators which properly correlate the positions of the electrons on neighboring chains 
without disrupting the singlet overlap property or the charge degrees of freedom greatly 
then the scaling dimension of pair hopping between the chains will be renormalized and pair 
hopping will become relevant and lead to an instability similar to that envisioned in the 
interlayer tunneling mechanism for superconductivity. Relevant operators having this effect 
occur naturally when the spin-spin superexchange interaction between chains is considered 
in the bosonization language [|19] . 
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FIGURES 

FIG. 1. Calculated Overlap for the ISE and Heisenberg Models 
Shown are the calculated overlaps between the N spin groundstate of the 1/r^, nearest 
neighbor Heisenberg and XY models and the states obtained by inserting a nearest neighbor 
singlet pair of spins into the N — 2 spin groundstates of those models. The dashed line 
is a fit to the ISE results using 0.0817 + 0.778A^~^, the sohd to Heisenberg results using 
0.820 + 0.740A^~^ and the dotted and dashed a fit to the XY results using 0.808 + .819A^-l 



FIG. 2. The Decay with Separation of the Overlap 
The exact overlap for the XY and the Monte Carlo calculated overlap for the ISE model 
as functions of the system size for spins inserted in a singlet configuration separated by half 
the size of the system. The results obtained in this way should be purely power law in the 
large system size limit, whereas correlation functions at fixed system size generally deviate 
from power law when the separation is not small compared to the system size. We have also 
examined correlation functions at fixed system size, and exact results for smaller systems for 
the ISE and Heisenberg models. All of those produced results consistent with a~2 behavior 
for the ISE and Heisenberg models and a~s behavior for the XY model. 
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